KYUSHU-HET-76 



00 



> 

o 

(N 
O 
(N 

O 



An asymptotic formula for models with caustics 



o 
o 

(N 

O ! Tomohiko Sakaguchi* 

Q '. Department of Physics, Kyushu University, Fukuoka 812-8581, Japan 



February 7, 2008 



Abstract 



We introduce an asymptotic formula for calculating quantum mechanical and quantum theoretical 
■ models with caustics, like the Nambu-Jona-Lasinio(NJL) model. This asymptotic formula is given 
by the form attached the extra term, which suppresses the divergence induced because of caustics, 
I to the leading term of the WKB approximation. This formula guarantees validity of the mean field 
^ ■ approximation of models with caustics. 
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1 Introduction 



There is the WKB approximation method as one of the nonperturbative approximation methods: 
the action is expanded around the classical solution, and terms with derivatives of greater than third 
order are treated as perturbations, to obtain a series expansion expressed in terms of a loop expansion 
parameter, such as h or I.Q. Then it is very important that positivity for all eigenvalues of a 
second derivative of the action at the classical solution, that is, the convexity condition of the action 
must be satisfied. When this approximation method is utilized to the Nambu-Jona-Lasinio(NJL) 
model, which is one of the most famous field theoretical models exhibiting dynamical chiral symmetry 
breaking phenomena, many interesting results are made in a level of the mean field approximation. 
For example, T. Inagaki et al. [2j have investigated the phase structure of the NJL model at finite 
temperature and chemical potential in an external magnetic field and gained the interesting results 
such as the magnetic catalysis. 

Nevertheless the mean field approximation in the NJL model is open to question: for example, in 
reference we had calculated the contribution of one loop diagrams of the auxiliary fields for the 
the effective potential and Gap equation of the NJL model. Then the infrared divergence occurs from 
the loop diagram of the massless auxiliary field. This fact implies that we always have to introduce 
a current mass in order to suppress this infrared divergence, and either do the mean field approxima- 
tion. Otherwise the WKB approximation fails because the convexity condition is not satisfied. This 
discourages the interesting results obtained in the chiral limit. Therefore it is important to investigate 
whether these mean field approximations is justified in terms of the asymptotic expansion. 

Generally when one and more eigenvalues of the second derivative of action vanish at the classical 
solution, higher order terms of the WKB approximation diverge. This situation is often called 'caus- 
tics'. L. S. Shulman 4 says that this divergence occurs approximately because the third derivative of 
action does not vanish at the classical solution. In fact, he has given the ad hoc treatment of caustics 
for a Green function in quantum mechanics. Unfortunately, This method cannot apply to calculate 
the effective potential for the NJL model in the chiral limit. We need to consider other method for 
the treatment of caustics. 

It is well-known that similar circumstance like caustics occurs when we calculate an asymptotic 
expansion of an integral using the saddle point method [S] [Ol- For example, let us consider the Bessel 
function 



where Re a stands for the real part of a. If a = 1, the second derivative of the exponent becomes zero 
at the saddle point z = 0. In such a case, there are two a-dependent and nearly coincident saddle 
points in the neighborhood of a = 1. Utilizing the two saddle points and replacing the exponent 
asinhz — z with a cubic function with respect to a new variable, C. Chester et al. have presented the 
method giving a uniform asymptotic expansion as a — > 1 and N ^ oo [7j [S] [S]. The application for 
an integral with multi- variables had been presented by N. Bleistein |lflj . 

We would like to apply their method to the NJL model. However this method cannot use directly 
to the NJL model, because the classical equation of the NJL model has only one classical solution or 
at least we do not know whether the another solution, which coincides with the classical solution in 
the chiral limit, exists. We need to make the method that we can use for the model having only one 
classical solution so that we can obtain the same property as the asymptotic expansion presented by 
C. Chester et al. 

In reference an other method obtaining an asymptotic expansion of an integral which the 
second derivative of the exponent of the integrand becomes zero at the saddle point is introduced. 
Although this method does not give a uniform asymptotic expansion, it needs only one expansion 
point and the result coincides with the one based on the method presented by C. Chester et al. in the 




(1) 



neighborhood of caustics. It is the most remarkable thing in this method that this expansion point 
always exists even if we do not know whether the another saddle point exists. Utilizing this fact, we 
introduce other approach for obtaining the asymptotic formula of the models with caustics. 

This paper is organized as follows: in section 2, we mention the method obtaining the asymptotic 
formula of an integral with one variable, which had been utilized in the paper This method 

utilizes a point at which the second derivative of the exponent becomes zero as an expansion point of 
the exponent, rather than the saddle point. Then we rewrite the obtained formula with respect to the 
saddle point and show the formula is the form attached the extra term which suppresses the divergence 
to the result of the saddle point method. In the end of section 2, we show that the formula coincides 
with the leading order of the asymptotic expansion given by C. Chester et al. in the neighborhood 
of caustics. In section 3, we apply this method to obtain an asymptotic formula of an integral with 
multi- variables. In the same way as the result of section 2, We show that the formula is the form 
attached the extra term which suppresses the divergence to the result of the WKB approximation. 
This implies that the mean field approximation of this formula coincides with the one of the WKB 
approximation and gives the important result that the mean field approximation of the model with 
caustics, like the NJL model, is justified in the sense of the formula whose result converges. The final 
section is devoted to conclusions and discussions. 



2 An asymptotic formula of an integral with one variable 

In this section, we establish an asymptotic formula as ^ oo of an integral 

I{a,N)= [ g{z)e^f^''"Uz, (2) 



where the real functions g{z), f{z,a) are analytic functions of their arguments, is a large positive 
parameter and c is some contour in the complex z plane. The function f{z, a) has the saddle point 
zq which is dependent on a. The second derivative on the steepest descent, d'^f{z,a)/dz'^\^^^g^^-^, 
becomes zero at some value of a, say a, and takes negative value for a > a} 

We expand f{z,a), g{z) around z which satisfy with /(^)(z,a) = (and f^^\z,a) = 0): 

I{a,N) ~ / {g{2) + 0{z - z)) 



c 



X exp 



AT J f{z, a) + f^'\z, a){z - ~z) + 1/(3) (i, a){z - zf + 0{{z - zf) 



dz 



(3) 



where we utilized the abbreviations /("^(z, a) = d'^f{z,a)/dz^. We assume that f^^\z,a) / 0. Con- 
sidering only in the neighborhood of a, we can regard f^^\z,a) as very small quantity, especially, we 
set that \ f^^'^{z,a)\ ~ 0(1/A^ 3) for sufficiently large A^. So we can consider both terms /^^^ (5, a) (z — 5) 
and f^^\z, a){z — z)^ as same order terms, in the neighborhood of z on the z space. We retain terms, 
g{z), f{z,a), f^^\z,a){z — z) and f^^\z,a){z — z)^ as the principle part, and remaining terms with 
higher derivatives are treated as negligible quantities, to obtain 



I{a,N) ~ gi~z) ( / dtea*--^*, (4) 



^In the case of the integral with one variable, it is expected that the another saddle point, say Zq, which coincides with 
the saddle point zq at a = a always exists. In this paper, as we would like to derive the method utilized only one saddle 
point, we assume that we cannot deform the integration path so that it traces the steepest descent of the point Zg when 
a > a and we do not utilize the point Zq. Also, we do not consider about a < a because the saddle points zq and Zq become 
the complex conjugate each other and it is out of our aim. 



where we have set 



a) 



"^'nV if ('Hz, a) 



/«(5,a)^C(«) 



and C' = iV^C- The integral part of Eq. just corresponds to the Airy integral 

1 



Ai(C') 



2iTi 



dt e3 



Therefore we can obtain the asymptotic formula 



I{a,N) ~ (27ri)<7(i) 



3 pNf{z,a) 



— Ai(C') 



(5) 



(6) 



(7) 



J(3)(5,a)y Nfi 

We can rewrite the asymptotic formula (PT|) with the saddle point zq, not the point z. When a is 
in a neighborhood of d, a distance between z and zq, \z — zq\, becomes very small. So we can expand 
g{z), f{z,a) and its derivatives at the point z around the saddle point zq: 



g{z) = g{zo) H , 

f{z,a) = f{zo,a) + ^f(^\zo,a){~z-zof + ^_f'^'Hzo,a){2-zof + -- 

f('H~z,a) = f(^){zo,a){~z-zo) + \f('Hzo,a){z-zof + ---, 

= /(2)(5,a) = fi^)(zo,a) + f('Hzo,a){z-zo) + --- , 
f('\~z,a) = /(3)(^o, «) + ... . 

If we assume that /(2)(zo ,a) 0(1/N3 ) and that N is sufficiently large, Eq. (fTT|) gives 



z- Zq 



1 / ' 

iV3 



/(3)(zo,a) 

and Eqs. ©, 0, liTH]) and (|T^ are approximately given by 

g{z) PS g{zQ) , 

/(z,a) « /(zo,a) - -— 77577 TT^ 

3 {f(-i){zo,a)y 



2 /(■^)(zo,a) 



O 



2 ; ; 



(8) 

(9) 

(10) 

(11) 
(12) 

(13) 

(14) 
(15) 

(16) 
(17) 



respectively. Substituting Eqs. (fTl)) - H17() into Eq. (pil) . We obtain the asymptotic formula of the 
integral in terms of the saddle point zq: 



Iia,N) ~ (27ri)5(zo) 



/(3)(zo,« 



I ^7V/{^o,«)-k'^(") 



A^3 



-Ai(C') 



27r 



and C{a) is deformed as 



C(«) 



iV/(2)(zo,a 



2(/(3)(zo,a))2 



.(2V^i)C'^e-i«^Ai(0, 



(18) 



(19) 



The form of the asymptotic formula (jlBj) is the form attached the extra term to the leading order of 
the asymptotic expansion based on the saddle point method, so that the divergence due to the term 
/(^^(zo, a) at a = d is suppressed. Expanding the Airy function © by use of the saddle point method 
as follows 



Ai(C') 




63 



\ 1 

4 



(20) 



and substituting this expansion in the formula H18j) takes off the extra term and leaves the leading 
term of the asymptotic expansion of the integral ((2j by the saddle point method. We find that, as in 
Schulman's textbook Ai, the term f^'^^zo, ct) plays a important role in order to converge the result of 
the WKB approximation. 

Finally, we show that the asymptotic formula (jlBf) coincides with the uniform asymptotic expansion 
of the integral (jSJ based on the method presented by C. Chester et al. in the neighborhood of a. The 
leading order of the asymptotic expansion due to their method is given as 

f,NA{a) 



where 



Aia) 
C(a) 

Go (a) 



I{a,N) ~ (27ri)ao(a)- 



-Ai(C') , 



2C5 



a) 



+ 5(^0)1 



20 



a) 



(21) 

(22) 
(23) 

(24) 



A point Zq represents the another saddle point that coincides with the saddle point zq when a ap- 



proaches to a? When a is in a neighborhood of a, a distance between Zq and zq, 



zq\, becomes 



very small. Therefore we can expand giz^), f{zQ,a) and its derivatives at the another saddle point 
Zq around the saddle point zq: 



9(4) 


= 9{zo) + 


+ i/(2)(zo,a)(zo*- 


^o)' + ^/=^)(^o,a)(^o-^o)' + --- 


(25) 




= f{zo,a) 


, (26) 


= f^'Hz*o,a) 


= f^'Kzo, 


«)(^o - ^0) + \f^'^ 


[zo,a){zQ - zof H , 


(27) 




= f^'Hzo, 


a)+/(3)(zo,a)(zo* - 


- zo)^ , 


(28) 


f^'H4,c^) 


= f^'\zo, 


a) + • • • . 




(29) 



If we assume that f^'^\zo,a) ~ 0(1/A^3) and that N is sufficiently large, Eq. (P7|l gives 



zo 



2/(2) ( 



zo,a) 



O 



1 



and Eqs. (P5|l and (|^ are approximately given by 



f{zl,a) ^ f{zo,a) 



2(-/(2)(zo,a))= 



3 (/(3)(zo,a))2 



f'\z*„a) 
f^'\4 



a) 



-f^'Hzo,a) 
f^^\zo,a) , 



(30) 



(31) 

(32) 
(33) 



^In making change of variables, the saddle points zq and Zq are corresponded to — ^2 and ^2 , respectively. 



respectively. Substituting Eqs. (ElJl-ljSSl) into Eqs (HU, and we obtain 

"l(-/(2)(zo,a))=^ 



A{a) w /(^o,a) - ^TTTsTT a\2~ ' 



2 (/(3)(zo,a))^ 



(34) 



ao(«) « ^^(-o) , (35) 

and substitution of these expressions in the asymptotic expansion (|21j) gives the asymptotic for- 
mula H18p . This fact shows that, even if we do not know whether the another saddle point exists, we 
can derive the asymptotic formula (|18|) consistent with the result of C. Chester et al. at least in the 
neighborhood of a. 

3 An asymptotic formula of an integral with multi- variables 

Next we investigate an asymptotic formula of an integral with n component variables, which is formed 
as 

/„(a,Ar) = [ e^^(^'")d"2; , (36) 



D 

where x = {xi, X2,- " > , d"x = dxidx2 ■ ■ ■ dxn and D is a region in n dimensional space. A saddle 
point F{x,a) is expressed as 

xo (a) = (xoi (a) , a;o2 (a) , • • • , xon (a) ) . (37) 
Without loss of generality an expansion of F{x, a) around the saddle point xo can be given by 



1 " 

F{x,a) = F{xo,a) + -'^Xi{a){xi - xoif 

1=1 

1 

+— ^ aijkia){xi - xoi){xj - xoj){xk - xok) H , (38) 

i,j,k=l 

where we assume that Ai(a) becomes zero at a = d and takes negative value for a > a, Aj(a) < (i = 
2, • • • , n) for a > a and am (a) ^ 0. We would like to give an asymptotic formula of the integral H36() 
by extension of the method discussed in Section 2. The most direct method is to find a point x so 
that F(x,a) is expanded as 



1 " - 

F{x,a) = F{x,a) + ai{a){xi - xi) + -'^Xi{a){xi - Xif 

i=2 

1 " 

+— ^ dijk{a){xi - Xi){xj - Xj){xk - %) H • (39) 



If we can find out such a point x, we obtain the asymptotic formula 



,.(„, ~ (2„)e-<-. (I) ' (^) ' (n {-^')) A.(C') . (40) 



where 



and we have only to rewrite this formula by F{xq, a) and its derivatives at the saddle point xq. 

Although such a point x always exacts, it is difficult to exactly find out the point. Fortunately, we 
do not need to derive explicitly the point. As the same way in section 2, when a is in a neighborhood 
of a, a distance between x and xq, \x — xq\, becomes very small. Therefore we can expand F{x,a) 
and its derivatives at the point x around the saddle point xq: 



F{x, a) 
Fi'\x,a) 

= FP{x,a) 

= Fil\x,a) 
Fi\x,a) 

FS(i,a) 



1 1 

F{xo, a) + -Ai(xi - xoif + — aiii(xi - XQif H , 

Ai(xi - xoi) + - xoi)^ H , 

Xi{xi - xoi) + ^aiii(xi - xoi)^ H , (i = 2, • • • , n) 



Ai + aiii(xi - xoi) + 

amixi - xoi) H = 

Xi6ij + aiij{xi - xoi) + 

am H , 



{i = 2,--- ,n) , 
= 2,--- ,n) , 



where we have used the abbreviations 



d^F{x,a) 



(42) 
(43) 

(44) 

(45) 
(46) 
(47) 
(48) 

(49) 



We assume that Ai(a) ~ 0{1/N'i) and that is sufficiently large. As in an integral with one variable, 
Eqs. dm), ^ give 



Ai 



Xi - Xoi 



am 



Xi Xoi 



am 



,n 



and Eqs. (|17|) and are approximately given by 

1 (-Ai)=^ 



F{x, a) 
Fm{x,a.) 



F{xo,a) - 
l(-Ai)^ 



3 al^^ 



2 am 
\i{a)5ij , (i,j 
aiii(a) , 



2,- 



(50) 



(51) 

(52) 

(53) 
(54) 



respectively. Although the terms F\^^ {x,a) {i = 2, ■ ■ ■ ,n) remain finite, these turns out to be sup- 
pressed as higher order terms. The use of Eqs. (|51|) - (|54j) gives 



In{a,N) 



dTx exp 



NlF{xQ,a) 



-Ai)^ 



3 af^i 



1 (-Ai)^ , _ . am , ~ \3 , 1 \ / 

O (Xi -Xl) + — — (Xl -Xl) +-\i{Xi 

I am 3! 2 



Xif + • • • 



Making changes of variables 



Xi — Xl 



r ( I ^1 ) -^i — ) — 2, • • • , TLj , 

A^3 V«iii/ A^2 



C(«) 



(-Ai(«)) 
2aiii(a) 



3x1 



(55) 

(56) 
(57) 



we obtain an asymptotic formula of the integral ()36() 



In{a,N) 



^NF{xo,a) 




n 

27r\ 2 



1 



(58) 




coincides with the saddle point a;o at a = a. Also, in the same way as section 2, it is not difficult to 
show that the asymptotic formula (|58|) coincides with the uniform asymptotic expansion derived by 
N. Bleistein ^Uj in the neighborhood of a. 

Finally, we can simply comment about the mean field approximation. As same way as the asymp- 
totic formula (|18|) , the form of the asymptotic formula (|58|) is the form attached the extra term to the 
leading order of the asymptotic expansion based on the WKB approximation. Therefore the mean 
field approximation based on the asymptotic formula (|58j) coincides with the one based on the WKB 
approximation. This fact gives the important result that the mean field approximation of the model 
with caustics is justified in the sense of the formula ((SHI)- 



In this paper, we derive the asymptotic formula for the integral with caustics, which we can utilize 
when the integral has only one saddle point xq or at least when we do not know whether the another 
saddle point, which coincides with zq at a = a, exists. This asymptotic formula coincides with the 
uniform asymptotic expansion in the neighborhood of d. Of course, if we can find two saddle points, 
we should utilize the method presented by C. Chester et al. [HIHlini and N. Bleistein ^0] so that we 
can obtain the uniform asymptotic expansion. Also, we show the important result that the mean field 
approximation of the model with caustics is justified in the sense of the asymptotic formula whose 
result converges. This result guarantees the validity for the mean field approximation of the NJL 
model. 

Our next aim is to apply this result to quantum field theoretical models with caustics. We think 
that carrying out this is very important: in reference [Hj, we calculated the effective potential and 
Gap equation of the NJL model. Then in order to suppress the infrared divergence which occurs from 
the loop diagrams of the massless Nambu-Goldstone boson, we introduced a fermion mass in advance. 
According to reference TP, the result of the WKB approximation becomes bad not only at a = a 
but also in the neighborhood of a, whereas the result of the approach mentioned in section 2 gives 
very accurate approximation in this region. Therefore if the introduced mass is sufficiently small in 
comparison with the typical energy scale, there is possibility that the asymptotic formula derived in 
this paper gives better result than the WKB approximation. In this sense, re-estimate for quantum 
(= loop) effects of auxiliary fields in the effective potential and gap equation of the NJL model is 
very interesting. Also, in reference ^2j, A. Osipov et al. first calculated the quantum correction of 
the gap equation and effective potential of the auxiliary fields for the NJL model with the 't Hooft 
interactions and dealt with the change of the vacuum state due to it. Then they suggested that 
including the quantum correction causes caustics at some values of the coupling constants. The use 
of the asymptotic formula to these caustics is interesting issue. 



4 Conclusion and Discussion 



Acknowledgements 

He would like to thank Koji Harada for reading this manuscript and making useful suggestions. 



References 



[1] T. Kashiwa, Y. Ohnuki and M. Suzuki, Path Integral Methods (Clarendon Press, Oxford, 1997). 

[2] T. Inagaki, D. Kimura and T. Murata, Prog. Theor. Phys. Suppl. Vol.153, 321 (2004), Prog. 
Theor. Phys. Vol.111, 371 (2004), ^^^0307289, 

[3] T. Kashiwa and T. Sakaguchi, Phys. Rev. D68, 065002 (2003), lhe p-th/ 0306008| 

[4] L. S. Schulman, Techniques and Applications of Path Integration (Wiley-Interscience, New York, 
1981). 

[5] E. T. Copson, Asymptotic expansions (Cambridge university press, 1965). 
[6] F. W. .J. Olver, Asymptotics and Special Functions (A K Peters, Ltd, 1997). 
[7] C. Chester, B. Friedman, and F. Ursell, Proc. Camb. Phil. Soc. 53, No. 3, 599 (1957). 
[8] B. Friedman, SIAM Journal 7, 280 (1959) 
[9] F. Ursell, Proc. Camb. Phil. Soc. 61, 113 (1965). 
[10] N. Bleistein, J. Math, and Mech. 17, No.6, 533 (1967). 

[11] T. Kashiwa and T. Sakaguchi, Prog. Theor. Phys. Vol.110, No.3, 589 (2003). 
[12] A. A. Osipov and B. Hiller, Eur. Phys. J. C35, 223 (2004). 



